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Express Mail NO. EL319729134US 

METHODS AND APPARATUS FOR HELICAL 
RECONSTRUCTION FOR MULTISLICE CT SCAN 

BACKGROUND OF THE INVENTION 

This invention relates generally to methods and apparatus for computed 
tomographic (CT) imaging, and more specifically to methods and apparatus for 
acquiring and reconstructing helically scanned, medical CT images using a multi-slice 
CT imaging system. 

5 In at least one known computed tomography (CT) imaging system 

configuration, an x-ray source projects a fan-shaped beam which is collimated to lie 
□ within an X-Y plane of a Cartesian coordinate system and generally referred to as the 

p= "imaging plane". The x-ray beam passes through the object being imaged, such as a 

f ! patient. The beam, after being attenuated by the object, impinges upon an array of 

10 radiation detectors. The intensity of the attenuated beam radiation received at the 
detector array is dependent upon the attenuation of the x-ray beam by the object. Each 
;L, detector element of the array produces a separate electrical signal that is a 

M measurement of the beam attenuation at the detector location. The attenuation 

> measurements from all the detectors are acquired separately to produce a transmission 

D15 profile. 

In known third generation CT systems, the x-ray source and the 
detector array are rotated with a gantry within the imaging plane and around the object 
to be imaged so that the angle at which the x-ray beam intersects the object constantly 
changes. A group of x-ray attenuation measurements, i.e., projection data, from the 

20 detector array at one gantry angle is referred to as a "view". A "scan" of the object 
comprises a set of views made at different gantry angles, or view angles, during one 
revolution of the x-ray source and detector. In an axial scan, the projection data is 
processed to construct an image that corresponds to a two dimensional slice taken 
through the object. In a helical scan, a table on which the object is resting moves so 

25 that the object itself moves though the imaging plane while it is being scanned. A 
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multi-slice CT imaging system has a plurality of parallel detector rows configured to 
acquire attenuation measurements corresponding to one or more two-dimensional 
image slices of an object. The number of image slices and the thicknesses represented 
by the slices is dependent upon how (and whether) attenuation measurements from the 
5 parallel detector rows are combined. 

One method for reconstructing an image fi'om a set of projection data is 
referred to in the art as the filtered back projection technique. This process converts 
the attenuation measurements from a scan into integers called "CT numbers" or 
"Hounsfield units", which are used to control the brightness of a corresponding pixel 
10 on a cathode ray tube display. 

ri Helical reconstruction algorithms for multi-slice CT have been a focus 

for many studies. Li one known CT imaging system, a reconstruction algorithm is 
iH implemented for two special helical pitches: 3:1 and 6:1. This algorithm utilizes two 

i2 conjugate samples fi-om different detector rows to estimate projection samples at a 

5^ 15 reconstruction plane using linear interpolation. Although this method performs 
=; satisfactorily in many cases, it has a number of shortcomings. First, the sampling 

pattern is not always optimum because only two samples on either side of the plane of 
I y reconstruction are selected. For example, samples that are closer to the reconstruction 

Q plane but located on the same side of the plane will not be utilized. Second, a 3:1 

'=^20 helical pitch is non-optimal for projection sampling, because the first and last detector 
rows measure identical ray paths, reducing the amount of non-redundant information 
that is acquired. In fact, in one known helical reconstruction implementation, 
measured projections (after calibration) of these two rows are summed first before 
reconstruction takes place. Third, sharp structures in the original object (along a z- 
25 axis) are suppressed and degraded slice sensitivity profiles are obtained because linear 
interpolation suppresses high fi*equency information in the sampled data. 

It would therefore be desirable to provide methods and apparatus for 
helical reconstruction in multi-slice CT imaging systems that overcome the above- 
described shortcomings of known image reconstruction systems. 
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BRIEF SUMMARY OF THE INVENTION 

There is therefore provided, in one embodiment of the present 
invention, a method for imaging an object with a computed tomographic (CT) 
imaging system that includes steps of heUcally scanning the object v^ith a multi-slice 
CT imaging system to acquire attenuation measurements of the object, the 
5 measurements including more than two conjugate samples for estimation of a 
projection at a plane of reconstruction of the object; and filtering and backprojecting 
the attenuation measurements of the object, including the more than two conjugate 
samples, to reconstruct at least one image slice of the object. 

This embodiment and others provide, among other advantages, an 
10 improved sampling pattern and better use of the attenuation samples obtained during a 
scan. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Figure 1 is a pictorial view of a CT imaging system. 

Figure 2 is a block schematic diagram of the system illustrated in 

Figure 1 . 

15 Figure 3 is an illustration of a sampling pattern for image 

reconstruction in an embodiment of the invention employing a 2.5:1 helical pitch 
scan. 

Figure 4 is an illustration of interpolation, extrapolation, and weighted 
interpolation-extrapolation techniques to estimate an attenuation value at a point x. 

20 Figure 5 is a graph of slice sensitivity profile measurements using a 

thin slice phantom, comparing a slice profile using a known technique utilizing a 3:1 
helical pitch scan to an embodiment of the present invention utilizing a 2.5:1 helical 
pitch scan. 
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DETAILED DESCRIPTION OF THE INVENTION 

Referring to Figures 1 and 2, a computed tomograph (CT) imaging 
system 10 is shown as including a gantry 12 representative of a "third generation" CT 
scanner. Gantry 12 has an x-ray source 14 that projects a beam of x-rays 16 toward a 
detector array 18 on the opposite side of gantry 12. Detector array 18 is formed by 
5 detector elements 20 which together sense the projected x-rays that pass through an 
object 22, for example a medical patient. Detector array 18 may be fabricated in a 
single slice or multi-slice configuration. Each detector element 20 produces an 
electrical signal that represents the intensity of an impinging x-ray beam and hence the 
attenuation of the beam as it passes through patient 22. During a scan to acquire x-ray 
10 projection data, gantry 12 and the components mounted thereon rotate about a center 
of rotation 24. 

Rotation of gantry 12 and the operation of x-ray source 14 are 
govemed by a control mechanism 26 of CT system 10. Control mechanism 26 
includes an x-ray controller 28 that provides power and timing signals to x-ray source 

15 14 and a gantry motor controller 30 that controls the rotational speed and position of 
gantry 12. A data acquisition system (DAS) 32 in control mechanism 26 samples 
analog data from detector elements 20 and converts the data to digital signals for 
subsequent processing. An image reconstructor 34 receives sampled and digitized x- 
ray data from DAS 32 and performs high speed image reconstruction. The 

20 reconstructed image is applied as an input to a computer 36 which stores the image in 
a mass storage device 38. 

Computer 36 also receives commands and scanning parameters from 
an operator via console 40 that has a keyboard. An associated cathode ray tube 
display 42 allows the operator to observe the reconstructed image and other data from 
25 computer 36. The operator supplied commands and parameters are used by computer 
36 to provide control signals and information to DAS 32, x-ray controller 28 and 
gantry motor controller 30. In addition, computer 36 operates a table motor controller 
44 which controls a motorized table 46 to position patient 22 in gantry 12. 
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Particularly, table 46 moves portions of patient 22 through gantry opening 48. During 
a helical scan, this motion occurs while scanning is taking place. 

hi one embodiment of the present invention, various shortcomings of 
known CT image reconstruction systems are overcome by using more than two 
5 conjugate samples for estimation of a projection at a plane of reconstruction (POR). 
This embodiment uses as many samples as the sampling pattem supports. The 
samples are located within a predetermined distance from the POR. Another feature 
of this embodiment is that non-integer pitch helical scans are performed for pitches 
that are numerically less than the number of detector rows, hi other words, if the 

10 number of detector rows is N, a pitch P\\ is used, where P is not an integer, and P<N. 
For example, a 2.5:1 pitch is used in one embodiment having four detector rows so 
that no duplicated samples are acquired. Yet another feature of this embodiment is 
the use of nonlinear interpolation techniques to preserve high frequency image 
components. Examples of suitable nonlinear interpolation include Lagrange 

15 interpolation and weighted interpolation-extrapolation, among others. Interpolation 
takes place by weighting projections prior to the filtered backprojection. 

hi one embodiment and referring to Figure 3, a system 10 having four 
detector rows Rl, R2, R3, and R4 is used with a 2.5:1 helical pitch. Sampling 
patterns are depicted as a function of projection view angle, detector angle, and 

20 detector row. y^l, y02,y53, and /3A represent view angles at which corresponding 
detector rows Rl, R2, R3, and R4 intersect a POR (not shown). These angles are 
arbitrarily selected, as only a relative angular span between them is important for any 
selected helical pitch. For this embodiment, an angular distance between adjacent 
rows that intersect the POR is 0.8;z: . In Figure 3, regions REGl, REG6, and REGl 1 

25 are examples of conjugate regions, because samples in these regions differ in their 
view angles by either ;r or 2 tt . Each projection sample at the plane of reconstruction 
is estimated based on conjugate samples selected from three conjugate regions. 

A "region 16" that would be located above REGl 5 could be included, 
since that region would also fall within a predetermined distance from POR. 
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However, "region 16" is excluded in this embodiment to provide improved temporal 
resolution. For each set of conjugated regions, three regions are included in this 
embodiment. If a "region 16" were included, the conjugate regions would include 
four regions 1, 6, 11 and 16. To provide the best temporal resolution, this 
5 embodiment seeks to formulate a smallest data set, in terms of view angle span, for 
image reconstruction. Including "region 16" would increase the entire angular span, 
as well as increase computation time and reduce uniformity. 

Weights depend upon divisions of the region. To preserve symmetry 
and other properties, in one embodiment, projections are estimated along a curved 
10 plane written as >9/ = 2.8;r - P[ -2k - I3[ = 1 .Itt - y, and = OAtt - y , where 

>^2' >^3 5 ^4 represent view angles in a curved plane for corresponding detector 
rows Rl, R2, R3, and R4, respectively. Each curved plane is defined by boundary 
conditions given by regions cited in equations (1) to (4). In this notation, y represents 
the detector angle. Weights for third-order Lagrange interpolation for rows Rl, R2, 
15 R3, and R4, respectively, are written as: 




(1) 




(2) 
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In one embodiment of the present invention, weighted interpolation- 
extrapolation (WIE), i.e., a combining of weighted interpolation measurements with 
weighted extrapolated measurements, is used prior to the filtered backprojection. 
Figure 4 illustrates an example having three sampling points that are located on both 
sides of a point x where the interpolation is to take place. (A different case in which 
one sample is located at left and two samples at the right is treated similarly.) For one 
known helical reconstruction algorithm, only x2 and x3 are used for linear 
interpolation represented by line 50. Weights t2 and t3 used for x2 and x3, 
respectively, are written as: 



t2 = 



^3 = 



x3- X 
x3 — x2 

X- x2 
x3 — x2 



(5) 



15 



In one embodiment of the present invention, to use points xl and x2 to 
estimate x, extrapolation is used as represented by line 52. Weights for xl and x2 
(denoted by el and e2) are written as: 



-7- 



15-CT-5344 




x2 — X 




(6) 



In this notation, xl, x2, and x3 gre x-axis coordinates of locations of 
measured signals, and x is the location of the signal to be estimated. The weights used 
in linear interpolation are calculated based >on a relative distance of the two points to 
the interpolation location. Combined wights for points xl, x2, and x3 (denoted by 
ql, q2, and qS, respectively) are writtearas: 



q\ = {i-tl-tl)e\, 
q3 = {t2 +t")t3. 



t2. 



(7) 



In equation (7), a is a parameter that adjusts relative strength and 
contributions of the extrapolation. The terms ql, q2, and q3 represent weights for the 
10 three measured samples to estimate a sample in the POR. 

This interpolation scheme overcomes the shortcomings of linear 
interpolation and enables a better preservation of high frequency information contents. 
Weighting functions for rows Rl, R2, R3, and R4, respectively, for this interpolation 
scheme, represented by line 54, are written as: 
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where 0^ ^ P~P\ = /3 -OAtt -\- y • 

In equations (8)-(l 1), wl, w2, w3, and w4 are derived weights for four 
rows, based on conjugate regions shown in Figure 3 and in equation (7). For each 
sample to be estimated on POR, three conjugate samples are selected, based on Figure 
3, and equation (7) is applied to determine their weights. The three conjugate samples 
could, for example, come from any three of Rl, R2, R3, and R4. 

To demonstrate the advantages of the use of the above reconstruction 
weights expressed in equations (8)-(ll), a phantom study was performed. In this 
study, a thin plate placed parallel to an x-y plane inside a 20 cm poly phantom was 
scanned with 120kV/200mA/1.25mm/lsec at both 3:1 and 2.5:1 pitches. Images were 
reconstructed every 0.1 mm to ensure adequate samples to map out slice profiles. For 
the 3:1 helical pitch, a known reconstruction algorithm was used. For the 2.5:1 helical 
pitch, the weights expressed in equations (8)-(ll) were used. Profiles of the thin 
plates (and therefore, a system slice sensitivity profile) were calculated for both cases. 
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To avoid effects of statistical noise due to finite x-ray photon statistics, each point on 
a profile curve is an average over a 21 by 21 pixel region (in an x-y plane) centered on 
the thin plate. Figure 5 depicts a slice profile A for a 3:1 helical pitch image 
reconstructed using the known reconstruction algorithm and a slice profile B for a 
5 2.5:1 helical pitch image reconstructed using the weight functions expressed in 
equations (8)-(ll). A much narrower slice profile is obtained and the peak intensity 
for the profile is higher when equations (8)-(l 1) are used. 

To more fully evaluate a WIE reconstruction embodiment of the 
present invention, a noise analysis was performed. For comparison, images in z of 
10 WIE were smoothed by a filter kernel such that WIE with smoothing provides the 
same slice profile as the known algorithm without smoothing. A kernel that satisfies 
this condition has been found to be an 11-point kernel with coefficients [0.5, 1, 1, 
1, 1,0.5]. 

Image smoothing (in z) for the WIE images was also performed for a 
15 case in which the phantom was uniform in z. The standard deviation of the smoothed 
image was determined to be 5.23 HU. The standard deviation observed using the 
known algorithm was 6.63 HU. This result indicates that for the same slice sensitivity 
profile, a nearly 37% mA reduction can be achieved, if the noise is to be maintained at 
the same level. In images corresponding to the location where the thin plate is 
20 located, WIE-weighted images with smoothing show less image artifacts than images 
reconstructed using the known reconstruction algorithm. 

In another embodiment, an equivalent result to a non-linear 
interpolation is achieved by multiplying a set of projection date by a set of weights. 
The weights, for a third order Lagrange interpolafion, are written as equations (l)-(4), 
25 while weights for WIE are written as equations (8)-(l 1). 

It will thus be seen that embodiments of the present invention provide 
better sampling patterns, acquire less redundant information, result in less suppression 
of sharp structures of object and better slice sensitivity profiles than known helical CT 
image reconstruction methods and systems. 
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Although the exemplary embodiments described herein and having test 
results presented herein are 4-slice embodiments, other embodiments of the invention 
are applicable to multi-slice CT imaging systems providing different numbers of 
image slices. For example, in other embodiments, 8-slice, 16-slice, 32-slice, etc. CT 
5 imaging systems are used. In addition, the invention is not limited to embodiments 
providing a 2.5:1 helical pitch. This pitch was selected for purposes of illustration 
only. Other embodiments of the invention are applicable to CT imaging systems 
providing different helical pitches, hi addition, the CT system described herein is a 
"third generation" system in which both the x-ray source and detector rotate with the 
10 gantry. Many other CT systems including "fourth generation" systems wherein the 
detector is a full-ring stationary detector and only the x-ray source rotates with the 
gantry, may be used if individual detector elements are corrected to provide 
substantially uniform responses to a given x-ray beam. 

While the invention has been described in terms of various specific 
15 embodiments, those skilled in the art will recognize that the invention can be practiced 
with modification within the spirit and scope of the claims. Accordingly, the spirit 
and scope of the invention are to be limited only by the terms of the appended claims. 
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